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Abstract 

The critical behavior of two-dimensional (2D) anisotropic systems with weak 
quenched disorder described by the so-called generalized Ashkin- Teller model 
(GATM) is studied. In the critical region this model is shown to be described 
by a multifermion field theory similar to the Gross-Neveu model with a few 
independent quartic coupling constants. Renormalization group calculations 
are used to obtain the temperature dependence near the critical point of some 
thermodynamic quantities and the large distance behavior of the two-spin 
correlation function. The equation of state at criticality is also obtained in 
this framework. We find that random models described by the GATM belong 
to the same universality class as that of the two-dimensional Ising model. The 
critical exponent u of the correlation length for the 3- and 4-state random- 
bond Potts models is also calculated in a 3-loop approximation. We show that 
this exponent is given by an apparently convergent series in e = c — ^ (with c 
the central charge of the Potts model) and that the numerical values of are 
very close to that of the 2D Ising model. This work therefore supports the 
conjecture (valid only approximately for the 3- and 4-state Potts models) of 
a super universality for the 2D disordered models with discrete symmetries. 

PACS numbers: 05.50.-hq, 05.70.Jk, 75.10.Hk, 75.40.Cx 
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I. INTRODUCTION 



The critical properties of two-dimensional random spin systems have been extensively 
studied in the last few years Two-dimensional (2D) systems are particularly interesting 

due to a variety of reasons. Firstly, there are numerous examples of layered crystals under- 
going continuous antiferromagnetic and structural phase transitions @,|^. More recently, 
2D and quasi-2D crystals have begun to be fabricated and studied thanks to advances in 
deposition techniques, with an enormous increase in the variety of physical phenomena to 
be investigated Perfect crystals, however, are the exception rather than the rule, with 
quenched disorder always existing in different degrees. Even weak disorder may drastically 
affect the critical behavior, according to the celebrated Harris criterion 0. Secondly, the 
conventional field-theoretic renormalization group (RG) approach based on the standard 
0^ theory in (4 — e)-dimensions, and as applied to study properties of disordered systems 
by Harris and Lubensky and Khmelnitskii does not work in 2D due to the hard 
restriction e<^l. Similar considerations apply to the (2 + e) low-temperature RG approach. 
Thirdly, from a theoretical point of view a most challenging problem is to establish the 
relationship between random models and the corresponding conformal field theory (CFT) 
describing these at criticality. 

Some early exact results concerning the 2D random-bond Ising model (IM) with a special 
type of disorder (where only the vertical bonds are allowed to acquire random values, while 
the horizontal bond couplings are fixed) have been obtained by McCoy and Wu |]TU[. This 



type of ID quenched disorder without frustration was shown to smooth out the logarithmic 
singularity of the specific heat; the frustrated case was considered by Shankar and Murthy 
||TT| . Many years ago Dotsenko and Dotsenko initiated some considerable progress in the 
study of 2D random bond IMs by exploiting the remarkable equivalence between this problem 
and the N=0 Gross-Neveu model. For weak dilution the new temperature dependence of 
the specific heat was found to become C^lnlnr, r = '^-^ being the reduced deviation 
from the critical temperature Tc. However, their results concerning the two-spin correlation 
function at the critical point were later reconsidered by Shalaev Shankar and 



Ludwig [0]. By using the RG approach as well as the bosonization technique these authors 
showed that the large-distance behavior of this function at criticality was the very same as 
in the pure case. Some convincing arguments in favor of the critical behavior of the 2D 
IM with impurities as governed by the pure IM fixed point had been given earlier by Jug 
15| . Recently, a good number of papers devoted to Monte-Carlo simulations of the critical 



behavior of the random Ising model have been published [|T6|. Most Monte-Carlo data are in 



good agreement with analytical results obtained in [|12]-|T^. It should be mentioned however, 
that these analytical results have been obtained by employing the replica method. This, 
on the one hand, is known to give reliable results only in the framework of perturbation 
theory. On the other hand, the mathematical legitimacy of the replica trick has not yet 
been established. Moreover, replicas (though being very useful and convenient) appear not 
to capture the essentials of nonperturbative effects in the close vicinity of the phase transition 
point (Griffiths phase) (see, for instance, Jl^)- The study of nonperturbative effects in the 
critical properties of random systems is however beyond the scope of the present paper. 

Here, we point out that there is some scope to extend the previous analysis for the 2D 
random IM to other discrete-symmetry systems. A very interesting problem consists in 
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considering minimal CFT models with c < 1 as perturbed by randomness. These models 
comprise the 3- and 4-state Potts systems as particular cases, and these have interesting 
applications to real 2D crystals 0. Because the critical exponent a is positive for all these 
models the critical behavior is governed by a random fixed point in agreement with the Harris 
criterion. Some years ago Ludwig [0,11^ and Ludwig and Cardy ||20| made an attempt to 
calculate perturbatively the critical exponents of the random 3-state Potts model. Their 
approach was essentially based on the powerful CFT technique. More recently, Dotsenko, 
Pujol and Picco obtained the critical exponents for the dilute 3-state Potts model in a 
two-loop approximation by exploiting the Coulomb Gas representation for the correlation 
functions and a special kind of e-regularization, where e stands here for the difference between 
the pure system's central charge value and the conformal anomaly for the pure 2D IM (|). 
Dotsenko, Dotsenko, Picco and Pujol [^] have also found the new universality class of 
the critical behavior as corresponding to the broken replica symmetry proposed by Harris, 
Dotsenko, Stinchcombe and Sherrington ||23|| . 

Another interesting possibility is to study critical phenomena in 2D dilute anisotropic 
systems with many-component order parameters. The analysis of the critical behavior of 
such systems in (4— e)-dimensions has been developed in great detail many years ago ||2^, but 
cannot be directly applied to the 2D case. Therefore it would be interesting and important to 
consider studying these 2D models. This is the main goal of our paper. The key ingredient of 
our treatment is a fermionization trick first suggested by Shankar |2^ for the N-color Ashkin- 
Teller model (see also [^,^]). This method is quite general and may be extended to other 
systems. The initial Landau Hamiltonian as written in terms of scalar fields can be shown 
to map onto a multifermion field theory of the Gross-Neveu type with a few independent 
quartic couplings. This transformation can be done for Hamiltonians containing only even 
powers of each order parameter component, the fourth-order term being an invariant of the 
hypercubic symmetry group (this is the so-called generalized Ashkin- Teller model (GATM)). 

The work presented in this paper is organized as follows. In Section II we consider in 
brief the critical behavior of the weakly-disordered 2D Ising model with random bonds, this 
being the central theme of this research field. The tranfer matrix formalism is set up and 
the corresponding equations are written down. The computation of the two-spin correlation 
function for pure and random models at criticality is also reviewed. In Section HI we give a 
description of the fermionization trick allowing us to study the critical behavior of the pure 
N-color Ashkin- Teller model. In Section IV the critical properties of two interacting N- and 
M-color quenched disordered Ashkin- Teller models are studied. The RG method is used to 
obtain the exact temperature dependence of the correlation length, specific heat, susceptibil- 
ity and spontaneous magnetization near criticality, as well as the two-point spin-correlation 
function and the equation of state at the critical point. In Section V, exploiting the approach 
of Dotsenko, Picco and Pujol and of Ludwig we compute the critical exponent of the corre- 
lation length in a 3-loop approximation for the weakly-disordered minimal models of CFT, 
in particular for the 3- and 4-state Potts models with random bonds. We find that while 
for the GATM the introduction of disorder leads to critical behavior as characterised by the 
random-bond IM fixed point, for the minimal models of CFT this Ising behavior, conjectured 



by a number of authors recently for the 2D Potts models ||28|, is actually only approximate. 
The accuracy with which the Ising values of the exponents is observed, however, justifies 
the use of the term "IM superuniversality" for all these models, when disordered. Section 
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V contains a discussion and some concluding remarks. 



II. THE TWO-DIMENSIONAL ISING MODEL WITH RANDOM BONDS 

A. Transfer Matrix, Effective Action and RG for Thermodynamic Functions 

We begin with the classical Hamiltonian of the 2D Ising model with random bonds 
defined on a square lattice with periodic boundary conditions: 



TV 



H 



i=l,j=l 



(2.1) 



where i,j label sites of the square lattice, Sij = ±1 are spin variables, Ji{i,j) and J2{i,j) 
are horizontal and vertical independent random couplings having the same probability dis- 
tribution, which reads: 



P{x) = (1 - p)6{x -J)+ p6{x - f) 



(2.2) 



Also, p is the concentration of impurity bonds and both J and J' are assumed to be positive 
so that the Hamiltonian favors aligned spins. Notice that both antiferromagnetic couplings 
(creating frustration) and broken bonds (J' = 0) lead to ambiguities in the transfer matrix 
and must be excluded in the present treatment. Let us now consider the calculation of the 
partition function of the model under discussion 



Eexp( 



'T' 



(2.3) 



where H is defined in eq. ( p.l| ) and the sum runs over all 2^' possible spin configurations. The 
partition function is known to be represented as the trace of the product of the row-to-row 



transfer matrices Tj [2S-31 



N 



Tr\{Ti 



(2.4) 



i=l 



The Hermitian 2 x 2 matrix Tj rewritten in terms of spin variables reads [p9| -|31|P]: 



fi = exp I Ji(^,j>3(j)(T3(j + l)j exp (^^E ■^2(^,0^^1(0 



N 



(2.5) 



where cXa, a = 1,2,3 are Pauli spin matrices; here J2 and J| are related by the Kramers- 



Wannier duality relation [2S-31 



tanh(^) = exp(-^; 



(2.6) 



In eq. ( p.5|) we have set an irrelevant factor to unity. Since the non-averaged operator Tj in 
eq. (|2.7| ) is random, the representation in eq. ( |2.4[ ) is in fact inappropriate for computing the 
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partition function. In order to get a more convenient starting point for further calculations 
we apply the replica trick. We introduce n identical "replicas" of the original model labeled 
by the index a, a = 1, . . . , n and use the well-known identity for the averaged free energy: 



F = -Tin Z = -T lim -(Z" - 1) (2.7) 



Substituting eq.(|2.4|) into eg. (^.71) one obtains: 



n N 



F = -T hmlTr I] 11 - 1}- (2-8) 

In contrast to the case of random-site disorder, for the random-bond problem the two ma- 
trices T" and Tj with different row indices i ^ j depend on two different sets of random 
coupling constants and commute to each other for any a and (3. This allows us to average 
these two operators independently. After some algebra one arrives at 

^ = Trf ^ (2.9) 
where the tranfer matrix T of the 2D random-bond IM is given by 0]: 



a=l 

{N T ^ Jl n 

Y: log[(l - p) exp(- Y: crm)^3ij + l) + Pexp(- Y: ^?0>?(J + 1))] 
j=l ^ a=l a=l 



N 

X exp ^ X! log 



7^E<(i) +Pexp(— $:<(j)) 

a=l / -'■ a=l 



(2.10) 



Setting p to zero (or J = J') one is indeed led to the well known expression for the T-operator 
of the pure IM 



T iV J* ^ 

fpiM = exp{-Ea3(j>3(j + l)}exp{-5:ai(j)} (2.11) 
^ j=i ^ j=i 

The T-matrix is known to possess the Kramers- Wannier dual symmetry. In the language of 
spin variables this nonlocal mapping reads p]|,pi[]: 

n{k) = asik)as{k + 1) T2{k) = iai{k)a-i{k) n{k) = J] ^^iM (2-12) 

m<k 

where the operators Ta{k) satisfy the very same algebra as the Pauli spin matrices <Ja{n). 
It is easy to see that if p = 0, |, 1 the T-matrix given by eg. ( |2.1CI| ) is invariant under the 
dual transformation. The plausible assumption that there is a single critical point yields the 
eguation for the critical temperature T^: 

2J' J 

exp(-— )=tanh(-) (2.13) 
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Notice that the point p = | is not the percolation threshold, because the coupling constants 
J and J' are assumed to take nonzero values with the ferromagnetic sign. Writing T in the 
exponential form 

f = exp{-H) (2.14) 

one obtains the partition function in the following form: 

Z = Trexp{-NH) (2.15) 

where by definition H is just the logarithm of the transfer matrix T (the "quantum" Hamil- 
tonian). In the thermodynamic limit N ^ oo the free energy is proportional to the lowest 
eigenvalue of the "quantum" Hamiltonian H p]|,pT]]: 

F = ~T\nTrexp{-NH) ^ NTEo H\0 >= Eo\0 > (2.16) 

Here |0) is the ground state of H which is assumed to be non-degenerate. Actually this means 
that we assume T > Tc. From eq.( p.lO| ) and eq.( |2.14| ) it follows that H is not a simple local 
operator. A crucial simplification occurs by taking the y-continuum limit with ay ^ (the 
lattice spacing along the y-axis). In other words, after calculating the logarithmic derivative 
of T with respect to ay and setting ay to zero the "quantum" Hamiltonian takes on the 
following simple form (for details see 0): 

r/lnT ^ f " 

H = -r-|K=o) = - E K,a^{j)a^{j + 1) + K, ^ 

""y j=l { a=l 

+ ^:K(J>3 (j + 1))' + Kit (2.17) 

a=l J 

The higher-order terms in the spin operators are known to be irrelevant in the critical region, 
so that they can be dropped in eq.( p.l7|) . The replicated Hamiltonian, eq. (|2.17| ), may be 
converted into the fermionic one by means of the Jordan- Wigner transformation POBTI]: 

m—l 

c"(m) = a^(m) n<(j)Q" 
i=i 

m—l 

a± = -{a,± ia2), = H 11 0'), a = l,...,n (2.18) 

^ /3=lj=l 

where c°'{m) and c°^(m) are the standard annihilation and creation fermionic operators 
which satisfy the canonical anticommutation relations: 

{c"(m),c^+(n)} = rt„„ {c"(m),c'^(n)} = (2.19) 

After making different species anticommute, the Klein factors drop out of H. For each 



species it is convenient to introduce a two-component Hermitean Majorana spinor field 33 
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ipi{n) = [c°(n) exp(— iy) + c"'^ (n) exp{i'^)] 



^^{n) = -^[c"(n)exp(2y) + c°''(n)exp(-4)] (2.20) 



./2o^ 4 
with standard anticommutation rules: 

{^:{n),^^{m)} = —d^^d^mn c,b= 1,2. (2.21) 
ax 

where ax is the lattice spacing along the x-axis. Using eq. (|2.21|) and the following relations 

a°{n) = 2c°^(n)c°(n) - 1 



a 



3{n)a^{n + 1) = [c" (n) - c"(n)][c" (n + 1) + c"(n + 1)] (2.22) 



one can easily rewrite the Hamiltonian, eg. ( [^.17| ), in terms of Majorana fermionic fields. 
Now let us notice that in the vicinity of Tc the correlation length ^ goes to infinity and 
the system "forgets" the discrete nature of the lattice. For that reason we can simplify the 
Hamiltonian by taking the continuum limit a^. — > 0. Perfoming simple but cumbersome 
calculations we arrive at the 0(n)-symmetric Lagrangian of the Gross- Neveu model IH 



d x[illJad^a + molpa-ipa + Mo(^aV'a) ] (2.23) 

where 7^, = d = 7^(9^, ^ = 1,2,^) = ip'^'jo and 

T — T 

mo r-. Ki - K2 T = " uo ~ fsTs + Ki (2.24) 

Here mo,Mo are the bare mass of the fermions and their quartic coupling constant, respec- 
tively. Notice, that if p <^ 1, Mq ~ p. Provided p = | and T = Tc we have Uq ~ {J — J'Y- 

The RG calculations in the one-loop approximation are very simple. In fact, the 0(n)- 
symmetric Gross-Neveu model being infrared-free in the replica limit n ^ 0, the one-loop 
approximation truly captures the essentials of the critical behavior of the model under 
consideration. The one-loop RG equations and initial conditions are given by: 

du in — 2)u^ d\n.F il — n)u 

u{t = 0) = uo; F{t = 0) = 1 (2.25) 

where u is the dimensionless quartic coupling constant, /3(m) is the Gell-Mann-Low function, 
7^^(m) is the anomalous dimension of the composite operator ipip = e(x) (in fact, the energy 
density operator), t = In— , A = is an ultraviolet cutoff and a and m are the lattice 
spacing and renormalized mass, respectively. Here F is the following Green's function at 
zero external momenta: 

F = ^ = / d^xd^y(^(xU(yU(OUm (2.26) 
dr J 
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The solution of these equations gives the temperature dependence of the correlation length 
^ and specific heat C in the asymptotic region t — oo, n = fl]]: 

TT -i 

u = — : F T 2- 
2V 

T 

C ^ j dtF{tf ~ Inln^ (2.27) 

These results follow from the solution of the one-loop RG equations, eq. ( p.25| ), but in 
fact it is worth noticing that they are a direct consequence of a renormalization statement 



valid to all orders in perturbation theory. Consider a version of the field theory, eq. (|2.23| ), in 
which the quartic term is decoupled by the introduction of a scalar Hubbard-Stratonovich 
field 0: 

/■ - - 11- 
L = j d^x[tPa{id + mo)Va + ^<P^ + ^go(f)i^ai^a] (2.28) 

with gooCy/uQ. As a consequence of the functional version of the classical equation of motion 



34 



5L 1 

T- = + l^goi^ai^a = (2.29) 

the vertex parts F of the correlation functions G^^,^^'^^ = {tlja{x)tpb{y)\'^c{z)ipc{z)) and 
^ab — {i'a{x)ipb{y)4>{z)) are hnked by the relationship 

r^''^^ = -^?or^^;^^^ (2.30) 

where it has been indicated explicitly that the quadratic insertion refers to the Oi = ipa'^a 
operator. Imposing the renormalization conditions, eq. (|2.30|) leads to |^| 



go/g= ^1^= z/^^Zu (2.31) 

where g is the renormalized coupling constant, Z,^ the 0- field renormalization constant and 
Zij is the quadratic-insertion renormalization matrix for the operators {Oi} = {ipip, 0^}. 



(0), 

Tplp 



u 



Since for n = we have = 1, eq. ( p.31| ) leads to the exact result (3^^\u) = — 2m7 
between the Gell-Mann-Low and the anomalous dimension functions, implying 271 //52 = — 1 
for the coefficients of the lowest-order nonzero terms in the expansion of these functions in 
u (that is: P{u) = ^2^^ + " " l-ip^iu) = jiu -!-•••). Solving, for instance, the RG equation 
for the specific heat function leads to the remarkable Dotsenko and Dotsenko's result for the 
leading asymptotic behavior, when r— s>0 

C-l'^llnxl'^i/^^-lnlnr (2.32) 

by virtue of the above n = exact results. Similar considerations lead to the announced 
behavior of the correlation length, ^. 
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The main conclusion of this Section is that the critical behavior of the 2D random bond 
IM is governed by the pure Ising fixed point. It implies that all critical exponents of the 
weakly disordered system are the very same as for the pure model. Randomness gives rise 
to the self-interaction of the spinor field which leads to logarithmic corrections to power 
laws. In the special case V = \ duality imposes strong restrictions, in particular it gives the 
exact value of the critical temperature which is believed to be unique. At the critical 
point the original lattice model and its continuum version described by the Gross-Neveu 
model Lagrangian become massless, irrespective of the value of n. We conjecture that there 
are only two phases divided by the single critical point given by the self-duality equation, 
eq. (|2.13| ). It implies that under this assumption the Griffiths phase shrinks to zero. 



B. Two-spin Correlation Function at Criticality 

In order to complete the calculation of the temperature dependence of other thermo- 
dynamic quantities we have to compute the susceptibility and spontaneous magnetization 
near T^. For these calculations we need to find the large-distance asymptotic behavior of 
the two-spin correlation function at criticality. The most effective way for calculating differ- 
ent correlation functions for the 2D IM is to use bosonization. Below we shall give a brief 
description of this procedure, exploiting simple physical arguments. 

Before recalling the principles of the bosonization method, however, let us show how a 
straight formulation of the problem in terms of pure fermionic fields leads to some difficulties 
even in the case of the calculation of the pure Ising model correlation function's exponent 
Tj (= i) at criticality. As shown, e.g., by Samuel the two-spin correlation function 



can be expressed in the lattice formulation as the partition function of a defective lattice 
where along the line Tqr of bonds joining the two sites (0, 0) and (0, R) the "bond strengths" 
Aj,= tanh( 72/7") must be replaced by A~^. Namely: 

G,(i?) = (sooW = Af(exp{-(A,-A,-i) ^ y\,y^,+l}) (2.33) 



were the lattice (y-) Grassmann variables {ylj^yij} have been introduced |]35| , |36| . After 



suitable transformations, leading to the quadratic term of the effective Grassmann action in 
eq. (|2.23| ) without replicas, and in the continuum limit, the above eq. ( p.33| ) reads 



Gy{R,T,) = X^,{exptTo dy^{0,y)ij{0,y)) (2.34) 

with To = (A~^ — Ac)/2Ac = 1 at criticality (Ac = — 1 for the isotropic model). The 
two-component Grassmann (or Majorana) field is the same as in eq.(|2.23D and is given by 



ip = a ^{y'^ y). A possible strategy [|I|] is now to evaluate the R—yoo behavior of In Gy{R,Tc) 
through an expansion in powers of Tq. Use must be made of the propagator {x = x^--^^) 

S,{x - x') = (7/i(x)z^(x'))o = ^[i- iVfK{x - x') (2.35) 

where /a is some cutoff function. The typical term in the expansion for \nG{R) involves 
the multiple integral 
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rR 

l2n{R) = / dyidy2 ■ ■ ■ dy2nTr[So{yi - y2)So{y2 - y^) ■ ■ ■ So{y2n - Vi)] 
Jo 

= anR + bnlnR + --- (2.36) 
from which the R-^oo critical correlator could be evaluated through 



\nG{R) = -J2 , ^2n{R) +R\nX, 



n=l 



An 



(2.37) 



(the odd- valued power terms vanishing). Taking the (conjectural) point of view that all 
terms in R must cancel exactly, the evaluation of the In R terms can proceed by taking 
the choice (natural, but leading to some ambiguities) /a = 1 and evaluating every other 
y-integral exactly 



hniR) 



R 



dyidy2 ■ ■ ■ dy-. 



2n 



R 



dyidy2 ■ ■ ■ dyn x 



lo (yi - 2/2) (2/2 - 2/3) ■ • • iy2n -yi) Jo 
ln[(l - R/y2)/{l - R/yi)] ln[(l - R/ys)/{l - R/y2)] ■ ■ ■ ln[(l - R/yi)/{l - R/yn)] 



(yi - y2)iy2 - 2/3) ■ ■ • (yn - yi) 

After a straightforward but laborious reparametrization of the integral , we arrive at 
InGiR) = R\nX, - ^ ^ °^ ^ hniR) 



(2.38) 



n=l 



hniR) 



dzidz2 ■ ■ ■ dzn-i . ^^^^ — n ■ 



sinhX^i Zi/2 \ sinh2;j/2 7a-i xiR — x 



R 



Rdx 



29„ InRA 



(2.39) 



with the i?-dependence now neatly factorized out and the cutoff A~a ^ conveniently rein- 
stated. The coefficient 6n is evaluated through the Fourier representation 



z/2 



sinh(2;/2) 



dp 
2^ 



Fip)e' 



■tpz 



Fip) 



cosh np 



(2.40) 



leading to 6*^ = I-!^ dp[F ip)]"- . Finally, we get (dropping the i?-terms) 



00 n ,rps^2n 

In GiR) = - E ^ (^) In Rk = -7] In RA 



(2.41) 



where 



v = Yl 



n=l 



n \n 



1 r+00 00 1 

— / dpY^- 

27r j-00 n 



7^2 



cosh Tip 



(2.42) 



The last sum converges to a logarithm and the p- integral can be evaluated, provided |To| < 1. 
For To = 1, eq.( |2.42| ) leads to 77 = 1/4 []1|; however, the standard prescription calls for 
To = -\/2 + 1 and this leads to a divergence in the summation. Clearly, this is associated 
with the use of a uniform cutoff function /a = 1, but it must be stressed that to date no 
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further progress in evaluating the spin-spin correlator at criticality, using solely the fermionic 
formalism, can be reported. The situation is even more delicate when disorder is introduced, 
thus the method of the fermionic tail Tqji must be abandoned. 

Let us now begin discussing bosonization, with the following action: 

L = J £x{iipdip + [mo + T{x)]ipip} (2.43) 

where ip is a. Majorana spinor and r(x) is a random Gaussian field with the following 
probability distribution: 



P[r{x)] ~ exp{-— / d'x[r{x)Y} < r{x)r{y) >= Uo6{x - y) (2.44) 

ZUq J 

In fact, the action eg. ( |2.43| ) describes free fermions moving in the random potential t{x), 
which in our case is responsible for local fluctuations of the critical temperature Tc in the 
dilute ferromagnet. After applying the replica trick and averaging over "all" possible config- 
urations of t{x) one gets the very same Gross-Neveu Lagrangian as given by eq. (|2.23|) . The 
representation of the square of the spin-spin correlation function of the pure 2D IM, that is 

G{x - y) =< a{x)a{y) > (2.45) 

in terms of the path integral over the real bosonic field of quantum sine-Gordon model 
was found by Zuber and Itzykson |Q (see also |^) and reads: 



Gix-yf = Z-'-^ f D0sin(v^0(x))sin(v^0(y))exp{-5} 

S = l f d'x{{d,<pf + ^ cos(v^0)} 
2 J vra 

Z = J D(f)exp{-S} (2.46) 

At criticality, mo = 0, the path integral being Gaussian, the result of its evaluation is easily 
seen to be: 

G{x -y) ^ \x- y\'^ (2.47) 

The representation for the two-spin correlation function may be extended to the dilute 
system by replacing the bare mass mo ~ r with the random one mo + t(x) into eq. ( p.46| ). 
Of course, in the inhomogeneous case the non-averaged G{x,y), being sample dependent, 
depends on x and y separately. The averaged correlation function G{x — y) at the critical 
point may be computed (even without using the replica trick) in two stages: (i) firstly, the 
square root of G{x,yY is formally evaluated by means of expanding it in a power series in 
t{x); (ii) secondly, the resulting expression is integrated with respect to t{x) (for technical 
details of the calculations see |^,^). The conventional RG equation for the renormalized 
averaged correlation function reads: 

d d 

^^djl + ^^""^d^ + viu)}GR{p,u,fi) = (2.48) 
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where /i is a renormalization momentum, l3{u) is the beta-function and r]{u) is defined as 
follows 

V{u)=m (2.49) 

The spin renormalization constant Zcr{u) and the renormalized correlation function are de- 
fined in the standard way: 



G{p, «o, A) = ZMGr{p, u, ii) (2.50) 

The Kramers- Wannier symmetry was shown to apply in some vanishing terms linear in u in 
the expansions for rj{u) and Z„{u) 0, that is: 

Z^iu) = I + 0{u^)- niu) =^- + 0{u^) (2.51) 

Given I3{u) and ri{u) in the one-loop approximation, the solution of the Ovsyannikov-Callan- 
Symanzik equation for the correlation function is quite simple: 

G{v) G{R) ~ R--^ (2.52) 

So, the Fisher critical exponent takes the very same value ?7 = | as in the pure model. Notice 
that in contrast to higher moments of the spin correlation function, the first one does not 
contain the logarithmic factor due to the above-mentioned dual symmetry [1^. From this 



remark it follows that the temperature dependence of the homogeneous susceptibility and 
spontaneous magnetization are described by power-law functions of the correlation length ^ 
(without logarithmic corrections like In^) 

M ~ ~ (-r)s[ln^^]T^ (2.53) 

The equation of state at the critical point may be obtained from the usual scaling relation: 

H ~ M— ~ M^^ (2.54) 
As we predicted, all critical exponents of the quenched disordered system are identical to 



those of the pure model, apart from some logarithmic corrections |15 



III. THE iV-COLOR ASHKIN-TELLER MODEL 



The N-color Ashkin- Teller model (ATM) was introduced by Grest and Widom and 
consists of a system of 2D Ising models coupled together like in the conventional 2-color 
model. The lattice Hamiltonian of the isotropic A^-color ATM reads: 
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N N 

H=Y.Hi{s'') + J, E E ^ae, 

a=l a^b=l <nn> 

N N 

= - E {Jj:^ts^ + MT.^>]?} (3-1) 

<nn> a=l a=l 

±1, a = 1,...,N, <> indicates that the summation is over all nearest- 
neighboring sites, Hi{s°') is the Hamiltonian of the pure 2D IM, ea = s^Sj is the density 
energy operator, and J4 is a coupling constant between the Ising planes. 

This model was shown to be the lattice version of a model with hypercubic anisotropy. 



describing a set of magnetic and structural phase transitions in variety of solids p4| , ^ . The 
corresponding Landau Hamiltonian reads: 



N N 

<f' = E'^a, (5.<f)' = E(W (3.2) 

a=l a=l 

where $ is an A^-component order parameter, ml ~ r, mq ~ J4 and Vq are some coupling 
constants. In particular, in the replica limit the Hamiltonian, eq. (|3.2| ), describes the random- 
bond IM (for Vq > 0,uo < 0). If fo = 0, a phase transition in the 0(A^)-symmetric 
model with nonzero value of the spontaneous magnetization is known to be forbidden by the 
Mermin- Wagner theorem If Vq^^O the spontaneous breakdown of the discrete hypercubic 
symmetry occurs at T,, > 0. Since the term with vq is strongly relevant, the perturbation 
theory expansion with respect to Vq is actually hopeless near Tg. 

By exploiting the operator product expansion (OPE) approach, Grest and Widom ob- 
tained the one-loop /5-function for the quartic coupling constant J4. If J4 < and N > 2 
the phase transition was shown to be continuous and the critical behavior belonging to the 



2D IM universality class 39 



The exact solution of the multi-color ATM in the large N-limit was found by Fradkin 
4^ , who developed a rather complicated formalism based on bosonic fields and showed that 
a second order phase transition with IM critical exponents occurs if J4 < 0. In fact, as 
was shown by Aharony |^^, the model with hypercubic anisotropy, eq. ( p.2|) , in the large 



A^-limit is equivalent to the IM with equilibrium impurities. Moreover, for the 2D case he 
predicted the Ising-type critical behavior with logarithmic corrections. Such being the case, 
one expects that the critical behavior is identical to the IM critical behavior. Since a = 0, 
Fisher's renormalization of the critical exponents [33 is inessential and gives rise only to 



logarithmic factors. Notice also that in contrast to the pure case the specific heat C is finite 
at Tc. The exact solution of the 2D IM with equilibrium defects obtained by Lushnikov HH 
many years ago confirms these conclusions. 

The effective method for solving the model under discussion, based on a mapping of 



the original model, eq.( |3.1| ), onto the 0(A^)-symmetric Gross- Neveu model, was suggested 
by Shankar ||2^ (see also [|26| , p7| ). In order to show this equivalence let us transform the 
partition function Z by applying the Hubbard- Stratonovich identity: 
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Z = J D$exp[-if($)] 

r /■ 1 1 1 ^ 1 

= J D^DXexp{-J d'x[-{d,^r + -ml^' + -t^o E + ^A(x)$2 + ^[A(x)]^]} 

= / DXexpi-^ [ d^xX^)[Zj[ml + tX{x)]f (3.3) 

where X{x) is an auxiliary field; Zj is the exact partition function of the 2D IM which is 
known to correspond to a path integral over Grassmann variables (Section II): 

Zi = j DipDipexp{- J d^xliipdip + kq^^]} (3.4) 



Now let us replace kq = mg in eq. (|3^) by ml + iX{x) and substitute eq.( ^^ into eq. (|3^ ). 
This replacement is based on the fact that the energy-density operator in the 0*^ theory is 
0^ while in 2D fermionic models this is given by ipil!. We have |25| , |27[| : 

N 



I DXexp{-^J d^xX^) j f{D^aD^a^ 



2uq 

exp{- / d'^xliipadipa + ("^o + iX{x))%l)ai'a]} 



/I 
l[DiJaDiJaexp{-SGN) (3.5) 

a=l 

where Sqn is the Gross-Neveu action, given by eq. ( ^.23| ). In going from eq. ( ^.2| ) to eq. (p. 51) 
it is assumed that uq has been rescaled in the following way: mq — >■ Mq = moo~^ so as to make 
Uo dimensionless (the prime will be ignored hereafter). We see that the discrete hypercubic 
symmetry of the N-color ATM evolves into the continuous 0{N) symmetry, hidden when 
the system approaches the critical point. 

The one-loop RG equations for the N-color ATM have been already obtained in Section 
II, these being eq. ( ^.25[ ) where we must set n = N. Solving these equations gives the 
temperature dependence of the correlation length and specific heat in the vicinity of the 
critical point p7| : 

1 r ,1,1 N-1 

e~r-i[ln(-)]: 



JV-2 



C~[ln(i)]5^ (3.6) 

T 

As for the calculation of the correlation function, one can apply the procedure described in 
Section II. Like for the random IM case, the term linear in u for i]{u) and Z„{;u) vanishes 
due to the Kramers- Wannier symmetry. This implies the anomalous dimension of the spin 
Sa to be equal to \. We get: 

G{R) ~ R--^ 

7 1 7(iV-l) 

X ~ r"4[ln-]4(iv-2) 
r 

1 1 (iV-l) 

M ~ (-r) 8 [In -] 8(^-2) 

H ~ M^^ (3.7) 
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Notice that these results are vahd only for > 2, J4 < 0. If J4 > 0, the discrete 75 symmetry 
ip — >■ 'j^ipjipip — i> —ipip is spontaneously broken. From the 75-symmetry breaking it follows 
that (ipip) 7^ 0. It means that we have a finite correlation length, or in other words, a 
first-order phase transition ||39| , |25[| . So, equations ( p.TD reproduce the well-known results for 



some particular cases: = 0, 1, cxd corresponding to the random-bond IM problem, Onsager 
problem, and IM with equilibrium impurities, respectively. 

The symmetric eight-vertex model (or Baxter model) is known to be isomorphic to the 
N = 2 color ATM in the vicinity of the critical line. The phase diagram of the 2-color ATM 
contains the ferromagnetic phase transition line beginning from the IM critical point and 
ending at the point corresponding to the 4-state Potts model. Along this line the model 
exhibits weakly-universal critical behavior, with the critical exponents continuously varying. 
For instance, the critical exponent a changes continuously from a = (IM) to a = | (4- 
state Potts model [^). The above results obviously show the special nature of the N = 2 
situation, due to the factor In this case the system under discussion is described by the 
0(2)-symmetric Gross-Neveu model, or, equivalently, by the massive Thirring model with 
the /^-function being equal to zero identically and presenting nonuniversal critical exponents 
p5| . Since the A^ = 3 color ATM is equivalent to the 0(3)-symmetric Gross-Neveu model 
which is known to be supersymmetric |^6|, this model should possess a hidden supersym- 



metry (see for details [p5| , p6| ). 

Notice that in contrast to the 2D case, the critical behavior of the A^-color ATM in 4 — e 
dimensions (0 < e < 2) is governed by either the Gaussian, or the cubic fixed point and 
never by the IM fixed point. The type of critical behavior crucially depends on the order 
parameter component number A^. If A^ > Nc{e), the RG flow arrives at the cubic fixed point; 
in the opposite case, A^ < Nc{e), the Heisenberg (isotropic) fixed point is stable. Here Nc{e) 
is the critical dimensionality of the order parameter, its expansion in powers of e being 



N^ie) = 4 - 2e - (^C(3) - ^)e' + 0{e') (3.8) 



where C(3) = 1.2020528 is the Riemann zeta function, Nc{l) = 2.9 [|8[. If e ^ 2,Nc 
decreases and all the cubic fixed points approach the IM fixed point, merging at e = 2, 
irrespectively of the value of A^ E7\ . 



IV. THE GENERALIZED ASHKIN-TELLER MODEL WITH RANDOMNESS 

Now we extend our study of the A^-color ATM to two interacting M- and A^-color 
quenched disordered Ashkin- Teller models, giving rise to a generalized Askin- Teller model 
(GATM). The Landau Hamiltonian is given by 



H 



^ n I 



ml + T2{x) 



2^2 



1 N 1 N+M 

^ a=l ^ c=N+l 



(4.1) 
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where ^k,k = 1, M + N is an (M + A^)-component order parameter, a = 1, ...,N,c = 
N+1, N+M, rriQ ~ r, f^, > Q and /i, z/ = 1, 2. Summation over indices in the quadratic 
operators is understood. We may study two types of impurities: (i) uncorrelated impurities, 
and (ii) correlated ones. In these cases the two-point correlators for the independent random 
Gaussian fields read: 

(^) {T^l{.x)r^{y)) = Zf,5^J{x - y), (ii) (r^(a;)r^(?/)) = Zo6{x - y) (4.2) 

In fact we study some multicritical point in the model under discussion, eq.( |4.1| ), since it 
has been assumed that mio = = fnQ ~ r. This model in 4 — e dimensions (without 
disorder) was initially studied by Bruce and Aharony and by Pokrovskii, Lyuksyutov 
and Khmelnitskii (without cubic anisotropy) [|50| and then in numerous other papers p4| . 



By applying the replica trick and the "fermionization" method described in Section III, one 
arrives at the following effective fermionic action involving several types of quartic fermionic 
interactions: 

H = J (fx {i^td'^t + ^O^fc^fc + + M2^?^c ^d^d + Wo'^aK'^c^c 

+ ziv&^-xj/^vp^f xj/f + + ro^'X^^c^^^} (4.3) 

where \E'^ is a (real) Majorana fermionic field, = l,...,n —>■ are replica indices, 

a,b = 1, ...,N, and c,d = N + 1, N + M. Naively one expects the appearance of two 
impurity quartic fermionic couplings in the replicated Hamiltonian, eq.( ^.3| ). Instead, we 
have one additional four-fermion vertex tq (absent in the bare action). This counterterm 
arises in the course of the renormalization procedure and provides the closedness of the 
operator algebra. In some sense the appearance of this term means violating the Harris 
criterion. The latter is indeed essentially based on the assumption of the existence of only one 
operator responsible for the impurity-induced interaction of the order parameter fluctuations. 

The one-loop RG equations for the six coupling constants Ufj,,Vu,r and w are given by 
(for n = 0): 

^ = -(N - 2)uj - 2ziui - Mw^ 

^ = _(M - 2)ul - 2Z2U2 - Nw^ 

^ = -w[{N - l)ui + (M - l)u2 + zi + Z2] 
dz 

— i = -2zi[zi + {N- l)ui + 2Mr] 

(JjV 

dz 

= -2z2[z2 + {M- 1)U2 + 2Nr] 

diT 

— = -r[{N - l)ui + (M - l)u2 + zi + Z2] - w[Nzi + MZ2] (4.4) 
The initial conditions for both (i) uncorrelated and (ii) correlated impurities are as follows: 
{i) ^i(O) = ^2(0) = ^0, r(0) = 0, {11) zi(0) = ^2(0) = 2r(0) = ^0 (4.5) 
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It is easy to see that if one sets M = N = 1 (random-bond Baxter or, equivalently, 2-color 



ATM) one arrives at the RG equations first obtained by Dotsenko and Dotsenko |51|. In this 
case the coupling constants decouple from the others; moreover, instead of two couplings 
we have only one coupling constant z = zi = Z2. It was shown that in almost all cases 
even weak disorder would drastically change the critical behavior of the 2-color ATM from a 
nonuniversal behavior to the Ising-type one, modified by some logarithmic corrections (see 
also P). Even though the critical exponent a of the pure model is negative for wq < 0, 
for uncorrelated defects we find that the critical behavior of the model under discussion is 
changed by the emergence of the new scaling field r. In the case of correlated defects with 



Wq < the critical behavior of the random model was shown |51|| to be still nonuniversal 
with critical exponents a and u depending on both wq and on the concentration of impurities 
(ro). This is the only exceptional case in which we would have nonuniversal critical behavior 
for a disordered system. In all cases the two-spin correlation function was shown to have, 
however, the same large distance behavior as for the 2D IM 

Now let us consider two interacting A^- and M-color ATM without randomness {z^ = 
To = 0). There are several different types of asymptotic behavior of the coupling constants 
Ufj,{t),w{t), but there is only one stable solution exhibiting infrared-free behavior. That is 
given by: 

111 

uJt) = — -, U2(t) = — -, w(t) = 0(- — ), t^oo (4.6) 

{N-2)t' ^' (M-2)t' ^' 4lnt^' ^ ' 

As a result, the original model decouples into two independent A^- and M- color models as 
described in Section III. Thus, the hidden symmetry of the model near the critical point is 
the continuous 0{N) x 0(M) group. There also exists a solution of the RG equations given 

by 

u^it) = U2(t) = ±-w{t) = — 1- -, t^oo (4.7) 

^' ^' 2 ^' {N + M-2)t' ^ ' 

corresponding to the higher symmetry 0{M + A^) being explicitly broken in the origi- 



nal Landau Hamiltonian, eq.( [4.1|) . This is shown to be unstable. For instance, provided 
A^ = 2, M = 1 (or viceversa) and 'Ui(O) = ^2(0) = \w{^), we would obtain the super- 
symmetric asymptotic solution of eq.( ^.7| ). Were these conditions to be broken, i.e. were 



the supersymmetry explicitly broken, this would not be restored in the infrared limit p6 
Notice that our model without cubic anisotropy and randomness was shown to exhibit this 
enhanced asymptotic symmetry in (4 — e)-dimensional space, provided M + A^ < 4 [Q. 

One may expect that, due to the critical decoupling of two multi-color ATM into two 
independent models, the quenched disorder does not affect the critical behavior of the system, 
eq. (|4.1|) . This is because if A^, M > 2 the specific heat is finite at criticality (eq.(^]5|)) and 
randomness is irrelevant in accordance with the Harris criterion. As was explained above 
this reasonable assumption should be checked in view of the obvious breakdown of the Harris 
criterion due to the appearance of the additional scaling field r. The answer is that this is 
indeed the case. In fact, it is easy to check that the solution given by eq.( [4.6|) and describing 
pure models is stable despite the presence of three disorder couplings. 

Thus, from our RG calculations it follows that, in contrast to a 2D IM with random 
bonds, weak quenched disorder here is irrelevant near T^. Moreover, in the critical region 
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the decoupling of two interacting multi-color ATM was found to occur even in the presence of 
quenched disorder. The temperature dependence of the main thermodynamic quantities near 
the critical point, the two-spin correlation function and equation of the state at criticality 
of the model under consideration are given by eq. (|3.5| ) and eq.( p.7|) . 



V. THE WEAKLY-DISORDERED MINIMAL CONFORMAL FIELD THEORY 

MODELS 

The critical behavior of the minimal models of conformal field theory with c < 1 and 
as perturbed by a small amount of impurities is far from being solved and therefore is 
of considerable interest. In accordance with the Harris criterion, weak quenched disorder is 
expected to be strongly relevant near criticality since the critical exponent a of these models 
is always positive and given by a = , with m = 3, 4, . . .. In particular, for the 3- and 

4-state Potts model we have a = |, (m = 5), and a = |, {ra = oo), respectively. 

The first results in this field were obtained in some pioneering papers by Ludwig 



by and Dotsenko, Picco and Pujol [^. They succeeded in developing a powerful approach 



closely connected with the formalism exploited in the previous Sections for describing the 
multi-color ATM. These authors suggested a special kind of e-expansion for computing the 
critical exponents, where now e = c — |. Here c is the central charge of the minimal models 
without randomness and \ is the conformal anomaly of the 2D IM. The main result of 
their considerations is that the fi{u) and 7^^(m) functions coincide with the corresponding 
functions for the 0(A^)-symmetric Gross- Neveu model obtained in the framework of the 
minimal substraction scheme combined with dimensional regularization. The distinguishing 
feature of this scheme is that these functions do not depend on e except for the first term 
in the /^-function. Thus, there is a clear possibility to apply the results of multi-loop RG 
calculations for the Gross- Neveu model in order to compute the critical exponents of random 
minimal models. At the present time we have the 5- loop expressions for the /?(M)-function 
and anomalous-dimension functions 7^(u) and 7^^(m) of the fermionic field ■?/' and composite 
operator t/^?/^, respectively, and as obtained in Unfortunately, these expressions contain 
a few unknown coefficients in the 4- and 5-loop terms. As for the anomalous dimension of 



the spin variable ^^(u), this function was obtained in in a 3-loop approximation. Notice 
that according to the conformal field theory classification, the spin variable corresponds to 
the operator $m,m-i5 whilst $1^2 = is the energy density operator. Thus, one may 
use the RG functions obtained only in the 3-loop approximation for the calculation of the 
critical exponents. 

Let us now compute the critical exponents of the correlation length and specific heat of 
the random minimal models in the 3-loop approximation. The corresponding expressions 
for the /9-function and temperature critical exponent function are given by ||52|| : 

^{u) = 2eu - 2{N - 2)u^ + A{N - 2)u^ + 2{N -2){N- 7)u^ 

7^^(m) = 2(A^ - l)u - 2{N - l)u^ - 2{N - 1){2N - 3)u^ 

3 — m , . 

e = — , m = 3,4,... 5.1 

2m 
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Here is the number of planes (colors), coupled to each other in the usual way like in the 
N-color ATM (eq. (|3.1| )). The critical behavior of the multi-color minimal models is governed 
by the nontrivial fixed point of eq.( |5.1| ). From this equation it follows that: 



11 1 e A^e^ 

- = - + = - + 2(iV - 1)1— + - ^^^] 

S[m — 1) 



where e takes on the discrete values defined in eq. (|5.1|) and z/q is the critical exponent of the 
correlation length of the pure model. 

To check the self-consistency of eq.( ^.l|) let us consider the limit A^ — oo, describing the 
system with equilibrium impurities. The result is easily seen to be: 

m + 3 1 — ao 

From the expression for the anomalous dimension of the order parameter r]{u) obtained in 



2T[ it follows that the Fisher critical exponent r]Q is unchanged in this limit, T]imp = t]q, 



where for the q-state Potts r/o is given by |53 



^0= , \ g = 4cos2^-, m = 2,3,5,cx) 5.4 

Am[m + 1) m + 1 

As was expected, we have obtained the duly renormalized critical exponent of the correlation 
length and an unchanged value of the order parameter anomalous dimension, in agreement 



with the predictions of the general theory |4^. As expected, in the trivial case A^ = 1 one 
gets the critical exponents of pure systems. 

The critical exponent values for random models can be easily obtained from eq. ( ^.2|) 
by expanding v^. in powers of e and setting A^ = 0. The results of these calculations are 
presented in Table I. The most striking feature of the above expansions is that they look like 
rapidly convergent series, even in the case of the 4-state Potts model. As a matter of fact, 
this is not so surprising, because in the Thirring model (A^ = 2) the anomalous dimensions 
7^^(m) and 7^(tt) are known to be some geometric progressions in u. Notice also that if 



we set A^ = 2 in eq.(|5.1|) we do not obtain these progressions The reasons why the 



0(2)-symmetric Gross- Neveu model in the minimal substraction scheme is not completely 
equivalent to the Thirring model have not been completely understood as yet . It is also 
important that all values of z/^ for any arbitrary integer m slightly exceed unity, so that all 
values of Ur are negative in full agreement with the Harris criterion. 

As was mentioned above, the 3- loop results concerning the critical exponent ri{u) were 



obtained in ||2T|. It turns out that the numerical values of rj are in the close vicinity of the 
Ising model value, rj = 0.25. The reasons for that are as follows. The 1-loop approximation 
term which is expected to give rise to the main contribution to the deviation of rj from rjQ 
vanishes due to the Kramers- Wannier dual symmetry. The 2-loop correction was shown to 
be proportional to ew^, therefore the deviations of r] from the pure values are proportional 

to e3 m. 
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Thus, the numerical values of critical exponents of the weakly-disordered minimal models 
are very close to the critical exponents of the 2D IM. This has clear implications for the 3- 
and 4-state Potts models with random bonds, as shown by Table I. From a numerical point 
of view one might be tempted to conclude that all these models are described by the 2D 
IM fixed point |2^. We see here that this "super universality" is only approximate, though 



accidentally verified to a high degree of accuracy. 



VI. CONCLUSIONS 

It has been shown that the critical behavior of a good number of 2D anisotropic systems 
controlled by the IM fixed point is stable in the presence of weak quenched disorder. This 
statement was found to hold quite generally for the 2D IM, multi-color ATM, and some of its 
generalizations for which randomness is marginally relevant. In the case of the 2-color ATM, 
or Baxter model, disorder drastically changes the nonuniversal critical behavior inherent in 
this model over to the Ising-type critical behavior. Although some of these models exhibit 
a breakdown of the Harris criterion, this does not affect, in general, the stability of the 
IM fixed point. It is commonly believed that the type of randomness (random-bond, or 
site-disorder) does not play a role near T^, despite the fact that random-site disorder has 
not been studied in great detail as yet. 

On the numerical side, Monte Carlo simulation results are in good agreement with the 



analytical results based on the RG calculations [|T^]. For instance, the high-accuracy MC 
simulation results for a 1024 x 1024 Ising lattice with ferromagnetic impurity bonds, recently 
obtained by Schur and Talapov [TB|, show that the exponent of the two-spin correlation 



function at criticality is numerically very close to that for the pure model. On the other 
hand, numerical results obtained by Domany and Wiseman |^ do somewhat contradict 



the theoretical predictions. These results, for a 256 x 256 lattice, favor a log-type behavior 
of the specific heat near Tc for the disordered 2-color ATM and 4-state Potts models, and 
a double-log behavior of the specific heat for the random-bond IM. As was established by 
Dotsenko and Dotsenko the specific heat of the random 2-color ATM should exhibit 



the double- log divergence at the critical point. 

The critical behavior of the random 4-state Potts model was shown to be described by 



a new fixed point which does not coincide with the IM one [pT] , p!9| . The conjecture made 
in that the perturbation theory expansion around the free fermion theory appropriate 
for the Ising model is valid till the end point of the ferromagnetic phase transition line 
(describing the 4-state Potts model) is actually incorrect. Exact results for the repulsion 
sector of the sine-Gordon theory obtained by means of the Bethe ansatz in |Q show that 



the perturbation theory around = An (free fermions) diverges at = ^ (see also | 
i.e. it has a finite radius of convergence, and cannot be continued to = Stt. 

The critical exponents corresponding to the disordered 4-state Potts fixed point slightly 
differ from the IM ones. The numerical results are believed to be sensitive to that difference. 
It is interesting to compare the estimate for the critical exponent u of the 4-state Potts model 
based on the 3-loop approximation with known numerical results. From Table I it follows 
that u = 1.081. Novotny and Landau [^] obtained for the Baxter- Wu model (equivalent 
to the 4-state Potts model) the following value: u = 1.00(7). The result of Andelman and 
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Berker is given by: v = 1.19. Finally, the recent result obtained by Schwenger, Budde, 
Voges, and Pfniir is as follows: u = 1.03(8). 

We end this section by giving a remark that logarithmic corrections to the power-law 
dependence and corrections to scaling may give rise to a dependence of effective critical 
exponents on the concentration of defects as observed in some numerical experiments 
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TABLES 



model 


m 


e 






1^2 






TIM 


4 


-0.125 


0.889 


0.099 


0.017 


0.003 


1.008 


3-PM 


5 


-0.2 


0.833 


0.139 


0.038 


0.008 


1.018 


TPM 


6 


-0.25 


0.8 


0.16 


0.052 


0.014 


1.026 


4-PM 


oo 


-0.5 


0.667 


0.222 


0.13 


0.062 


1.081 



TABLE I. Critical correlation length exponent for random minimal models: TIM (Tricritical 
Ising Model), 3-PM (3-statc Potts Model), TPM (Tricritical Potts Model), and 4-PM (4-state 
Potts Model), m denotes the minimal model, e = (3 — m)/2m, uq is the homogeneous exponent, 
fr = + vi + 1^2 + t's the random one and Vn denotes the n-loop contribution to Vr 
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